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Abstract 

The pion electromagnetic form factor F 7T (Q 2 ) is calculated at the leading 
order of pQCD using the running coupling constant (~ as(Q 2 {l — x)(l — y)) 
approach. The resummed expression for F n (Q 2 ) is found. It is shown that the 
effect of the infrared renormalons can be taken into account by scale-setting 
procedure as(Q 2 ) — ► as(e^® ^Q 2 ) in the leading order expression. 
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1. It is known that one of the important problems in the perturbative QCD 
(pQCD) calculations is a choice of the renormalization scheme and scale in a trun- 
cated perturbative series [1,2]. Of course, when working to all orders in as(Q 2 ) the 
choice of scale and scheme is irrelevant, because a measurable physical quantity fis 
the same for all choices. But this freedom can be a source of difficulties, when one 
truncates a perturbative series for f. Namely, a bad choice of the coupling constant 
scale may result in large higher-order corrections in such truncated series. The 
method proposed by Brodsky, Lepage and Mackenzie (BLM) in Refs.[l,2] for de- 
termination of the coupling constant scale seems to be a good tool for the handling 
of these problems. The BLM scale-setting procedure amounts to absorbing certain 
vacuum polarization corrections appearing at higher-order into the one-loop QCD 
running coupling constant. In other words, the essence of this method is to find 
an average virtuality of a gluon in a loop diagram and to use it as the scale in the 
running coupling constant. As was emphasized in [1,2], this method may be applied 
to both inclusive and exclusive processes. 

Recently, the BLM method has been generalized by performing the calculations 
with the running coupling constant as{— k 2 ) at the vertices of Feynman diagrams, 
where k is the momentum flowing through the virtual gluon line [3,4]. This approach 
is equivalent to a resummation of an infinite number of terms in the perturbative 
series for the quantity f. 

In our previous works [5,6] we have used the running coupling constant approach 
for calculation of the pion and kaon electromagnetic form factors at the leading 
order of pQCD. As a result we have found the perturbative series in as(Q 2 ) for 
Q 2 F M (Q 2 ). In this letter we shall resum the series obtained in Ref.[5] for the pion 
form factor using the technique of the Borel transform and show that an effect of 
such resummation can be taken into account by redefining the scale of as{Q 2 ) in 
the leading order expression. 

2. Let us start from the factorized QCD expression for the pion electromagnetic 
form factor 

FAQ 2 ) = t f dxdy<PAy,Q 2 )TH{x, y] Q 2 ,as{Q 2 ))<PAx,Q 2 ), (1) 
JO Jo 

where Q 2 = — q 2 is the four momentum square of the virtual photon 7*, Q 2 is 
the renormalization and factorization scale. In (1) terms of order m 2 /Q 2 , as well 
as higher Fock state components of the pion wave function have been neglected. 
Here the pion wave function (p w (x,Q 2 ) is non-perturbative quantity and gives the 
amplitude for finding the quark (antiquark) carrying the longitudinal fractional 
momentum x(l — x) within the pion. The function Th(x, y; Q 2 , cts{Q 2 )) is the hard 
scattering amplitude for the subprocess qq' + 7* — > qq' and can be found by using 
pQCD, 

T H (x, y- Q 2 , a s (Q 2 )) = a s {Q 2 )T l H {x, y, Q 2 ) \l + a s (Q 2 )T 2 H (x, y, Q 2 ) + ..]. (2) 
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The leading order term as(Q 2 )Tjj(x,y,Q 2 ) here is given by [7] 

as mn(x, y, 2 ) = g^f^y C F = 4/3. (3) 

The next-to-leading order correction T^(x,y,Q 2 ) was found in [8]. The one-loop 
QCD running coupling constant a s (Q 2 ) in (2), (3) is defined as 

a s (Q 2 ) = ^ , (4) 

where 

2 

Po = 11 - gW/, 

and 71/ being the number of quark flavors. 

In Ref . [1] the authors argued that in a meson form factor calculations the scale 
in (3) has to be chosen equal to the square of a momentum transfer of the exchanged 
gluon Q 2 = Q 2 {l — x){l — y). But it is evident that as(Q 2 (l — x)(l — y)) suffers from 
singularities associated with the behavior of a$ in the soft region x — > l,y — > 1. 
In other words, integrations over x,y in (1) extended to x — > l,y — > 1 indicate the 
appearance of non-perturbative effects. Therefore, until now in all pQCD calcula- 
tions the scale Q 2 was taken equal to some fixed quantity (as a rule, Q 2 = Q 2 /4) 
and results for F n (Q 2 ) were obtained in the so-called "frozen coupling constant" 
approximation. In our previous paper [5] for calculation of F n (Q 2 ) in the running 
coupling constant approach we have expressed as(Q 2 (l — x)(l — y)) by means of 
the renormalization group equation in terms of as(Q 2 ) and calculated obtained in- 
tegrals over x,y as principal value ones (see, [9,10]). For Q 2 F 7r (Q 2 ) we have found 
the following perturbative series in as(Q 2 ), 



\2 

n=l 



00 'asW 



Q 2 F^Q 2 ) = ^rU) 2 Y.[^ 1 ) S n , (5) 



where r and S n are the model dependent coefficients. In calculations we have used 



the asymptotic (f> a sy(x) and the two-humped Chernyak-Zhitnitsky <ficz( x , Q 2 ) wave 



functions [11], 

(x) 



'asy\ 




4>cz{x,Q 2 ) = <p asy {x){l + -C^(2x-l 



MQ 2 Y 


72 /Po | 







(6) 



Here f n = 0.093 GeV is the pion decay constant, C^ 2 (2x — 1) is the Gegenbauer 
polynomial, 72 = 50/9 is the anomalous dimension, /j,q = 0.5 GeV is the normal- 
ization point. The phenomenological pion wave function (f>cz(x,Q 2 ) depends on 
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the scale Q 2 . But in this article we ignore this dependence and take it equal to 
4>cz(x, Hq). The coefficients S n of the perturbative series (5) and r for (f) asy (x) an d 
<t>cz(x, /Xq) are given by the following expressions, respectively 

S n = Cn/%- 1 = (n - 1)! (n - 2 + (3^\ r = 1, (7) 

and 

r 14 / 25 64 1 
5„ = c„/5 n - 1 = (n-1)! +n 1 + ^ + — — + 

" a: v v ' 3 V 2 n+ 3 n + 4«— i 



25 8 35 
+ — + 



2 n-l 3 n 3 . 4n-l 



r = 5 (8) 



The ^ in both cases are proportional to (n — 1)!, as a result the series (5), 
(7), (8) are ill-defined and must be resummed using some procedures for removing 
divergences. The procedure for such resummation is well known; perform the Borel 
transform in (5) with respect to the inverse coupling constant (see, for example, 
[3,12]) and after that calculate the inverted integral applying some regularization 
methods. 

3. The Borel transform of the series (5) is 
B[Q 2 F](u) = Y,7—r V] c n . (9) 

71=1 \ n L >- 

Eq. (9) can be formally inverted to give 

\Q 2 FM)\ eS = (167 T /?r)2 T^exp \-Anu/f3 a s (Q 2 )} B[Q 2 F](u). (10) 

In these cases, when one finds integral (10), it defines the Borel sum [Q 2 F n (Q 2 )] res 
of the series (5). For the solution of our problem we may follow this general recipe 
to find [Q 2 F n (Q 2 )] res . But by means of the identity 

roc I 

Eq.(l) can be written in a more convenient way [5,6]. Indeed, after substituting the 
leading order term in the renormalization group equation's solution [10] 

aJ0 2 (l-x)(l-v)) asiQ2) 1-^1 (12) 

into (1), and changing of the variables x, y to z = ln(l — x), w = ln(l — y), 
introducing the formula (11) and integrating over z,w for the asymptotic wave 
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function we find 



1 

(2^ 



duexp — 4nu/ p as(Q 2 ) 



2 2 
+ ^ — 



1 - « 



it 



It is evident that 
S[g 2 F] as2/ (u) 



+ 



I 



2 2 
+ - — 



u 



(1-u) 2 (2-u) 2 1 
The Borel transform of the series (5), (8) is 

1 25 64 



u 



(13) 



(14) 



B[Q 2 F] cz (u) 



It is easy by using 



16 



14 



(2-u) 2 (3-u) 2 (A-u) 2 

50 8 140 
+ t; + 



3(1 -u) 2-u 3-u 3(4-«)" 



(15) 



d_ 

du 



n-l 



B[Q 2 F](u) 



u=0 



to recover the series (5), (7), (8). 

The Borel transforms (14), (15) have singularities on the real u axis (at u = 
1,2,3,4), which are known as infrared renormalons. The infrared renormalons are 
responsible for a factorial growth of the expansion coefficients in Eqs.(7),(8). This 
fact confirms our conclusion in [5] concerning an infrared renormalon nature of 
divergences in the perturbative series (5), (7), (8). 

The difficulties with the double and single poles at u — 1, 2 in (13) can be 
avoided by taking the integral in its principal value. Then we obtain 



Q 2 F n (Q 2 



(16tt/ t ) 



The similar calculations lead to 

(SOtt^) 2 



Q 2 F n (Q 2 



res 
CZ 



A) 



- + (In A - 2) -U- + (In A + 2) -j^ 



233 , U.Li(\) 



(16) 



25(ln A - 2) + 8(8 In A + 1) + 



A 2 

A/A , , 35 N L?(A 4 ) 
4(41nA + y)^^ 



A 3 



(17) 
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In (16), (17) Li(X) is the logarithmic integral [13] 



Li(\) = / A = Q 2 /A 2 . (18) 

jo m a; 

The principal value prescription used above to regulate the infrared renormalon 
singularities at u — 1, 2, 3, 4 in Eq.(lO) (or at x — 1 in (18 )) produces an ambiguity 
in the perturbative series for the pion form factor F n (Q 2 ). Any other approach for 
defining the integral (10) would differ in the treatment of these singularities. The 
ambiguity introduced by our treatment, as has been shown in [5], is a higher twist 
and for Q 2 F n (Q 2 ) behaves as A 2 /Q 2 . It is well known that appearance of infrared 
renormalon singularities in the perturbative calculations indicate an importance of 
non-perturbative higher twist effects [14]. A nature of the non-perturbative higher 
twist corrections is process dependent [9,10,14]. In the case of the pion form factor 
they, probably, may be associated with contributions from the pion's higher Fock 
state wave functions | qq'g >, etc. 

4. Some comments are in order concerning these results. First of all, it is in- 
structive to compare sources of the infrared renormalons in our example and in other 
QCD processes considered in [3,4]. In these articles the running coupling constant 
approach was used in one-loop order calculations for resummation any number of 
fermion bubble insertions in the gluon propagator. This technique corresponds to 
partial resummation of the perturbative series for a quantity under consideration. 
In fact, in pQCD an infrared-safe and renormalization scheme invariant quantity 
has the perturbative series in a s {^ 2 ) similar to that shown in (5). But now the 
coefficients S n of such series have the following expansion in powers of f3 

S n = c n ^- 1 + r n ^- 2 + - (19) 

The Borel transform of such perturbative series defined as in (9) (that is, using 
only c n from (19)) corresponds to partial resummation of the series (5) and the 
only source of terms of order ag/^ -1 in (5) is the running coupling constant [3,4].) 

We have emphasized in Sect. 2 that in our calculations we use only the leading 
order term as(Q 2 )T^. Therefore, it is not accidental that S n in (7), (8) have exactly 
~ Pq -1 dependence and, hence, the expressions (16), (17) are exact sums of the 
corresponding series. The source of term ~ agP^' 1 in (5) is the running coupling 
as(Q 2 (l — x)(l — y)), but here it runs due to the integration in (1) over the pion 
quark's (antiquark's) longitudinal fractional momentum. May in our case terms 
~ ag^Q^ 1 be also described by fermion vacuum polarization insertion into the 
exchanged between the pion constituents gluon propagator or not, require further 
investigations. 

5. It is useful to find how much the resummed expressions (16), (17) differ from 
the "bare" ones. In Figs.1,2, the ratio R wJ . = [Q^Q 2 )]^. / [Q 2 FAQ 2 )t 
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where [Q 2 F n (Q 2 )Y wf is the pion form factor in the frozen coupling approximation, is 
depicted. As is seen, the contributions of the infrared renormalons to [Q 2 F n (Q 2 )]™ S j 
are considerable for both wave functions. Nevertheless, in both cases these contri- 
butions can be hidden into the scale of as(Q 2 ) in as(Q 2 )Tj I (x,y,Q 2 ). Indeed, let 
us look for a function f(Q 2 ) — C\ + C20.S (Q 2 ), such that the leading order term 
asiQ^T^x.y.Q 2 ) with Q 2 = e fiQ2) Q 2 will give for Q 2 F 7T (Q 2 ) approximately the 
same results as Eqs.(16),(17). Numerical fitting allows us to get (rif = 3), 

f(Q 2 ) ~ -5.6 + 10.91« 5 (Q 2 ), for as „ 

f(Q 2 ) ~ -11.76 + 38.09a 5 (Q 2 ), for ^ C z- (20) 

It is worth noting that a dependence of f(Q 2 ) on a s (Q 2 ) was introduced and argued 
in the next-to-leading order scale-setting calculations in [2]. 

The results of numerical calculations of R w .f. with redefined scale in the leading 
order expression are shown in Figs.l, 2. As is seen from Figs.l, 2, R w .f. ~ 1 in the 
whole range of Q 2 . 

The redefinition of the scale Q 2 — > e^ Q ^Q 2 means, that an average virtuality 
of the exchanged gluon in the pion form factor diagrams equals to e^Q 2 . The 
coupling constant ag(e^ ^Q 2 ) for the functions under consideration (6) varies from 
a s ~ 0.39 at Q 2 = 2 GeV 2 to a s ~ 0.49 at Q 2 = 20 GeV 2 for CZ , and from 
a s ~ 0.54 at Q 2 = 2 GeV 2 to a 5 ~ 0.35 at Q 2 = 20 GeV 2 for asj/ . These values 
are typical for such kinds of studies [3,4]. 

It is important to notice that the functions (20) cannot be obtained directly 
from perturbative series (5), (7), (8) by means of the BLM next-to-leading order 
scale-setting prescription [2]. 

6. In this letter we have found the resummed expression for the pion electro- 
magnetic form factor. For regularization of the infrared renormalons induced diver- 
gences in the perturbative series for F n (Q 2 ) we have used the method of the Borel 
transform and the principal value prescription. We have demonstrated that the re- 
summed result depends on the pion model wave functions used in calculations. Our 
investigations of the infrared renormalon effects in the pion electromagnetic form 
factor calculations allow us to conclude that exclusive processes have two sources 
of infrared renormalon contributions; 

1) integration in an exclusive process amplitude over a hadron constituents' longi- 
tudinal momenta, 

2) the running coupling constant in higher-order Feynman diagrams. 

The self-consistent treatment of these two sources of the infrared renormalon 
effects, as well as their study by taking into account a dependence of (p^x ', Q 2 ) on 
Q 2 , are subject of forthcoming publications. 
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FUGURE CAPTIONS 



Fig.l. The dependence of the ratio R asy on Q 2 . Here the curves 1 and 2 correspond 
to R asy with [Q 2 F v {Q 2 )}° asy in the frozen coupling constant approximation and after 

scale-setting procedure Q 2 — > e^^Q 2 , respectively. In calculations the QCD scale 
parameter A has been taken equal to 0.1 GeV. 

Fig. 2. The same as in Fig.l, but for Rcz- 
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